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We introduce the concept of multi-point propagators between linear cosmic fields and their non- 
linear counterparts in the context of cosmological perturbation theory. Such functions express how a 
non-linearly evolved Fourier mode depends on the full ensemble of modes in the initial density field. 
We identify and resum the dominant diagrams in the large-fc limit, showing explicitly that multi- 
point propagators decay into the nonlinear regime at the same rate as the two-point propagator. 
These analytic results generalize the large-fc limit behavior of the two-point propagator to arbitrary 
order. We measure the three-point propagator as a function of triangle shape in numerical simu- 
lations and confirm the results of our high-fc resummation. We show that any n— point spectrum 
can be reconstructed from multi-point propagators, which leads to a physical connection between 
nonlinear corrections to the power spectrum at small scales and higher-order correlations at large 
scales. As a first application of these results, we calculate the reduced bispectrum at one-loop in 
renormalized perturbation theory and show that we can predict the decrease in its dependence on 
triangle shape at redshift zero, when standard perturbation theory is least successful. 
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I. INTRODUCTION 

Detailed understanding of the development of gravi- 
tational instability is one of the central issues for the 
study of structure formation in observational cosmology. 
Although the details of the formation of objects are af- 
fected at small scales by the presence of baryonic matter, 
in the context of a dark matter dominated universe it 
is thought that the global properties of the large-scale 
matter distribution is determined to a large extent by 
the collisionless dynamics of dark matter. Therefore, the 
Vlasov equation, i.e. the collisionless limit of the Boltz- 
mann equation, describes the dynamics we are interested 
in (see e.g. [J, [l3| for details). This is for instance what 
pure N-body cosmological simulations attempt to solve. 
At large scales, where orbit crossing is insignificant, the 
Vlasov equation reduces to the dynamics a pressureless 
perfect fluid. 

Although the problem is well posed, its resolution be- 
yond the linear regime is still largely an open question. 
Some ideas have been developed in the literature (see 
e.g. [l| and references therein) but to make predictions 
beyond the linear regime for the density power spectrum, 
the standard approach in the literature has been to resort 
to semi-analytic prescriptions obtained by fittin g to nu- 
merical simulations. These fitting formulae, e.g. [l6l.l23t. 
originate either from the near universal transform ad- 
vocated in Q or are based on an even more empirical 
construction, the halo model 



These prescriptions however offer predictions for the 
power spectrum with relatively low accuracy, at the level 
of 5-10%, and are insecure in cases of non standard cos- 
mological models. With the advent of precision obser- 
vations of the large-scale structure, in particular in the 
context of the search for dark energy signatures in the 
growth of structure, this is clearly insufficient [l3|- It 
then becomes crucial to build reliable, and well con- 
trolled, prescriptions for describing the nonlinear evolu- 
tion of the power spectrum and higher order correlation 
functions such as the bispectrum. For the latter there is 
not much published work on the expected behavior be- 
yond tree order, although fitting formulae analogous to 
the power spectrum case have been proposed 113] ■ 

Recently, a promising new approach to fitting measure- 
ments of statistics in N-body simulations and interpolat- 
ing to other cosmologies has been developed in 0]. 

With the Renormalized Perturbation Theory (here- 
after RPT) formalism introduced in a new approach 
emerged for the construction of perturbation schemes 
applicable to such observables. It has been success- 
fully applied to the two-point propagator Q, and to 
the density power spectrum in the weakly non-linear 
regime Q. These results show in particular how the 
acoustic peaks are affected by the mode-coupling effects 
in the growth of structure. This reformulation of the 
gravitationalpcrturbation theory also finds alternative 
approaches [H [11 El HE HI, which give further in- 
sights into resummation schemes introduced by RPT. 
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The aim of this article is to show that some of the core 
results that have been originally discovered for the two- 
point propagator in Q can also be obtained in the case 
of multi-point propagators, thus extending exact known 
results regarding the evolution of cosmic gravitational 
instability. 

The plan of the paper is the following. We first suc- 
cinctly recall the general formalism and the results that 
have already been obtained at the two-point level [3], 
involving the cross-correlation between fields at the ini- 
tial and final conditions. As shown in previous studies, 
the behavior of this correlator at large-fc is of crucial im- 
portance and provides us with solid ground for the ex- 
pansion of the correlation functions of interest @. We 
then show that the three-point propagator, namely the 
cross-correlation between the final density field and two 
of the initial modes, can also be computed explicitly in 
the high-fc limit. This quantity reflects how the final 
density is sensitive to the two mode-mode couplings. We 
then extend this result to propagators of arbitrary num- 
ber of points. This is the central result of this paper. 

To contrast our results with numerical simulations, we 
extend the method developed in Q to measure multi- 
point propagators. This allows us to compute the three- 
point propagator in numerical simulations against our 
predictions in the large-fc regime. We finally comment 
on the relevance of these results for the construction of 
power spectrum and bispectrum resummation schemes, 
and provide a first application to the case of the bispec- 
trum. 



II. DYNAMICS 

Cosmological gravitational instability in a pressureless 
fluid in the single stream limit (when orbit crossing can 
be neglected) is governed by a set of three coupled equa- 
tions relating the density contrast 5 = p/ p — \, the pecu- 
liar velocity field v and the gravitational potential. Con- 
sistently with this limit one can assume the fluid to be 
irrotational and thus fully described by 5 and the veloc- 
ity divergence, 6* = V • v. Following [l^, the equations 
of motion in Fourier space can be written in a compact 
form with the use of the two component quantity 



(<5(k,. 



(1) 



where the index a = 1,2 selects the density or velocity 
components, Ti. = dlna/dr is the conformal expansion 
rate, a(r) the cosmological scale factor and r conformal 
time. The time variable s is defined from the scale factor 
by 



the scales we are interested in to a very good approxima- 
tion by replacing a(T) by the linear growth factor of the 
cosmology under consideration [l| . 

The equations of motion then read (we henceforth use 
the convention that repeated Fourier arguments are in- 
tegrated over), 

9,^'„(k,s) + a,.b«'f,(k,s) = 

7ahc(k, ki, k2) ^'6(ki, s) ^'c(k2, s), (3) 



where 



-1 
-3/2 1/2 



(4) 



and the symmetrized vertex matrix ^abc describes the non 
linear interactions between different Fourier modes and 
is given by 



7222 (k, ki,k2) 
7i2i(k, ki,k2) 



(5D(k - ki - k2) 



fe(k-ki-k2) 



|ki+k2|2(ki-k2) 



(ki 



ki 



2fc2 



(5) 



7ii2(k, ki,k2) = 7i2i(k, k2,ki), and 7 = otherwise; 
with (5d denoting the Dirac delta distribution. The for- 
mal integral solution to Eq. ^ is given by (see [1, [l^, H^] 
for a detailed derivation), 

*a(k, s) = gafc(s) 0b(k) + / ds' gab{s-s') 



X7^^|,(k, ki, k2)vl/,(ki, s')*d(k2, s'), (6) 

where 0a (k) = (k,s = 0) denotes the initial condi- 
tions, set when the linear growth factor a(T) — 1 and 
s — Q. The linear propagator gab{s) is given by 



9ab[S) 



e° 



3 2 
3 2 



-3s/2 
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3 -3 



(7) 



for s > 0, whereas gab{s) for s < due to causality, 
and 5ab(s) ~* <5ab as s ^ 0+. We will work with growing 
mode initial conditions, as usually assumed in numerical 
simulations, for which 



t>a(M) = Ua<5o(k), 



(8) 



and Ua = (1, 1). 

A perturbative solution to Eq. ([6]) can be obtained ex- 
panding the fields in terms of the initial ones. 



fa(k,5) = £vi/H(k,5), 



(9) 



s — Ina(T), 



(2) 



and corresponds to the number of e-folds of expansion. 
We will first consider a cosmology with = 1 and 
f^A = 0. The extension to a general case is obtained at 



with 



*i")(k,s)= J d3ki...d3k„ <5D(k-ki...„) 



in) 

;6i bn...b, 



(ki, . . . , k„; s)0bj (ki) ...(t)b„ (k„ 



(10) 
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FIG. 1: Diagrammatic representation of the series expansion 
of (k) up to fourth order in the initial conditions (f). Time 
increases along each segment according to the arrow and each 
segment bears a factor gcd{sf — Si) if Si is the initial time 
and s/ is the final time. At each initial point and each vertex 
point there is a sum over the component indices; a sum over 
the incoming wave modes is also implicit and, finally, the time 
coordinate of the vertex points is integrated from s = to the 
final time s according to the time ordering of each diagram. 
For instance, at fourth order there are two different possible 
topologies. 

where are fully symmetric functions of the wave- 

vectors that can be obtained recursively in terms of gat 
and 'jabc [H • Note that these functions have a non-trivial 
time dependence because they also include sub-leading 
terms in e*. Their fastest growing contribution is of 
course given by the well known {Fn, G„} kernels in PT 
(assuming growing mode initial conditions), 

.Fi") = exp(ns) {^^„(ki, .., k„), G„(ki, .., k„)} (H) 

for a — 1,2 (density or velocity divergence fields respec- 
tively). 

The relation in Eq. ([6]) can be used to construct a dia- 
grammatic representation of the fields 4''^"). As an illus- 
tration, we show in Fig. [1] diagrams corresponding to the 
local field expansion up to fourth order. 

A detailed description of the procedure to draw the 
diagrams and compute their values can be found in Q, 
we can briefly summarize these rules here as follows. In 
Fig. [T] the open circles represent the initial conditions 
(/)b(k), where b — 1 {b — 2) corresponds to the density 
(velocity divergence) field, and the line emerging from it 
carries a wavenumber k. Lines are time-oriented (with 
time direction represented by an arrow) and have differ- 
ent indices at both ends, say a and b. Each line represents 
linear evolution described by the propagator gabisf — Si) 
from time Si to time s/. Each nonlinear interaction be- 
tween modes is represented by a vertex, which due to 



quadratic nonlinearities in the equations of motion is the 
convergence point of necessarily two incoming lines, with 
wavenumber say qi and q2, and one outgoing line with 
wavenumber q — qi -I- q2. Each vertex in a diagram 
then represents the matrix 7abc(qj qi, q2)- ft is further 
understood in Fig. [T] that internal indices are summed 
over and interaction times are integrated out in the full 
interval [0, s]. 

Note that at the level of the field \t', all diagrams are 
tree diagrams with no loops. This is unlike in standard 
field theories, and is due to the lack of "antiparticles" 
(both in the strict sense of the word, like in quantum field 
theory, and at the quasiparticle level, as e.g. holes in con- 
densed matter systems). However, loop diagrams appear 
once we calculate statistical averages such as correlators 
between fields. Other unusual properties compared to 
standard field theories are the non-trivial fc-dependence 
of vertices (see also fis'l) and last, but not least, a strong 
breaking of time-translation invariance. In fact, the the- 
ory is unstable with perturbations growing as a function 
of time from their initial conditions (represented by den- 
sity and velocity fields after decoupling). 

Having written the final conditions of interest 5* in 
terms of initial conditions (f), to build diagrammatic rep- 
resentations for statistical quantities such as the non- 
linear propagator or the power spectrum, one needs to 
specify the statistical properties of the initial conditions. 
Here we will assume Gaussian initial conditions, in which 
case the statistical properties are encoded in the initial 
spectrum of fluctuations, 

{M(l)Mn')) = SD{q + q')Pab{q). (12) 

with Pab = UaUhPo for growing-mode initial conditions, 
i.e. if Eq. ^ is satisfied. Ensemble averages are there- 
fore obtained by gluing together the open circles in the 
diagrams by pairs (according to the calculation rules for 
Gaussian fields) to form a symbol (g), which will appear in 
all diagrams in the figures below. This symbol represents 
the initial power spectrum Pab{q), when the outgoing 
lines away from the initial conditions have respectively 
indices a and b and carry wavenumbers q and q' . Finally, 
in all diagrams it is understood that internal wavenum- 
bers (corresponding to loops) when such symbols appear 
are integrated over. 



III. MULTI-POINT PROPAGATORS 

A. The Two-Point Propagator 

The idea at the heart of RPT is to realize that the di- 
agrammatic series described above for correlation func- 
tions can be partially resummed. This amounts to gen- 
eralizing the operator gab, encoding the linear evolution, 
to a fully nonlinear one, denoted by Gab, that effectively 
describes time evolution when coupling between modes is 
present. This generalized operator, called the nonlinear 
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propagator, is formally defined as, 



(13) 



This quantity basically measures how linear is the tran- 
sition from initial to final fields. Note that, although 
certainly 5'a(k) depends on (/'&(k') for k' 7^ k, transla- 
tion invariance demands that the expectation value of the 
derivative in Eq. is only non-zero if k = k' |23| . 

Due to nonlinear evolution, one expects Gab to decay 
to zero at small scales where non-linear effects become 
important and make ^'(k) very different from a linear 
transformation of 06 (k). In fact, for Gaussian initial con- 
ditions it is possible to give an alternative expression for 
Gab, as a cross-correlation between final VPa and initial 
(j)b fields i, 

(*a(k, s) M^')) - Gac{k, s) (0,(k) 0b(k')). (14) 

This equivalence gives rise to the alternative name of 
two-point propagator and reinforces the idea that decor- 
relation of fields '^a with respect to (pb mandates that Gab 
must vanish at small scales. Indeed, it was shown in Q 
that Gab can be computed in the large-A; limit and it fol- 
lows a Gaussian decay with k. This was also checked in 
numerical simulations (see Fig. [Hand discussion in sec. lIVI 
below). At low-A:, on the other hand, linear theory is re- 
covered and Gab(k ^ 0) = gab, as expected. 

Introducing the concept of a two-point propagator 
helps to understand nonlinear corrections to the power 
spectrum. The power spectrum can be written as Q 



P{k,T) = G^k,T) Po{k) + PMc{k,T), 



(15) 



where we have explicitly included the time dependence, 
Pq is the initial density power spectrum, and Pmc is 
the mode-coupling power spectrum. Note that what 
we call G in this equation is the density propagator 
G = Gil + G12 — GibUb, since we assume growing-mode 
initial conditions; a completely equivalent equation holds 
for velocities. The first term describes all the contribu- 
tions to P{k) that are proportional to the initial power 
spectrum at the same scale k, while the second, denotes 
the contribution to P{k) from modes other than k in 
the initial spectrum (hence called "mode-coupling"), i.e. 
Pmcik) is nonzero even if Pq vanishes at wavenumber k 
(as long as Pq is not zero everywhere). 

The first term in Eq. (jlSp has the most direct informa- 
tion of the initial power spectrum, since it is proportional 
to Pq at the same scale, thus varying k one can probe 
Po(fc) as long as the first term is dominant, i.e. at large 
scales where G has not yet decayed significantly. The 
two-point propagator G can be given an alternative in- 
terpretation by noting that Eq. p4)) says that {G/D+)'^, 
where I?+ is the linear growth factor, is the fraction of 
the power that grows by linear theory at a given k. As G 
decays, a given mode looses memory about its primordial 



(linear) value, and an important fraction of the power is 
contributed by nonlinear effects alone, i.e. Pmc- How- 
ever, that does not mean that information on the initial 
power spectrum is completely lost at such scales. In fact, 
RPT gives a precise expression for Pmc in terms of convo- 
lutions over the initial spectrum (see Fig. 5 in Q for how 
well this works when compared to simulations). How- 
ever the information on Pq contained in Pmc is rather 
crude, since PMc(fc) depends on nonlinear combinations 
of Pq (q) 's from a wide range of scales q weighted by com- 
plicated convolution kernels. 

As we will show in the rest of the paper, these kernels 
are the (square of) multi-point propagators and, in fact, 
they also determine the tree-level higher-order spectra 
(bispectrum, trispectrum, etc). Therefore, multi-point 
propagators will lead us to a different way of writing the 
mode-coupling contributions to P(fc) (and higher-order 
spectra), and a direct physical connection between non- 
linear corrections to the power spectrum at small scales 
and higher-order correlations at large scales. 



B. Definition of Multi-Point Propagators 

Having discussed the importance of the two-point 
propagator let us introduce the concept of multi-point 
propagators, denoted F'^^'^ as a natural extension of the 
definition in Eq. JTc 



for the three-point case and. 



<5^*a(k, S) 



p\ ^6(bbA^i)...5cj,bM 



) 



<5D(k-ki...p) Fi^,; (ki,...,kp,s),(17) 

where ki . p = ki -|- . . . -I- kp, for an arbitrary number of 
points. Note that F'^p^ denotes the {p + l)-point prop- 
agator, by translation invariance it depends only on p 
wavenumbers in Fourier space. 

The question we want to address in this section is the 
calculation of F*^^'^. We will shortly show that exact re- 
sults can indeed be obtained for F^^^ in the large k limit. 
This is a nontrivial generalization of the two-point case 
and constitutes a core result in this paper. In addition, as 
anticipated above, we will show in section lVlTl that corre- 
lation functions can be reconstructed solely in terms of F 
functions and initial power spectra Pq, e.g. for the power 
spectrum we will show that. 



P(fc) = E ^7 -^DCk - qi...r) [rW(qi, . . . , q.) 

xPo(gi)...Po(g.) d3qi...dV, (18) 

where F^^^ is the two-point propagator G used in Eq. (fT51) 
and the mode-coupling spectrum Pmc is described here 
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by the terms with r > 2. Equation (|18p is something 
already found in (J| within the Zel'dovich approximation, 
but extended here to the exact dynamics. 

Formally, the F functions can be computed order by or- 
der in a field expansion and thus they can be represented 
diagrammatically. Examples of such diagrams are shown 
in Fig. [51 To lowest order it can be easily shown from 
Eqs. (|9I10I17P that F^p) coincides with the kernels in 
Eq. ((T0| . Diagramatically this lowest order corresponds 
to a tree level diagram. Let us explicitly compute F^^^ 
starting from the second diagram in Fig. [T]and the rules 
described in Sec. |TT] to get an expression for "^^^^ that 
after the functional derivatives in Eq. pT]) gives, 

rifc!;trcc(ki,k2,s) = / ds' 5ad(s-s')7de/(ki,k2,k) 
Jo 

X5e6(s')g/c(s')- (19) 

For the density field (a = 1) and growing mode initial 
conditions this leads to, 

^(2) /, , ^ 9« /5 kix k2X 2x'^\ 

ri.Uc(ki,k2,.).... = ^"(7+2fc^+2fc7 + — j 

^ / 3 kix k2X 2x^ \ 



+ e 




with a; = ki • k2/fcifc2. The fastest growing term in 
Eq. (|5D)) is precisely the F2 kernel in standard FT. 

We now show that the equivalence in Eq. (fT4|) between 
two-point propagators and cross-correlations for Gaus- 
sian initial conditions, can be extended to multi-point 
propagators. Let us start with the three-point propaga- 
tor F(2). 

f 2") 

On the one hand, an expression for F)^^^ (ki, k2) can 
be formally written in terms of the functions JF^") using 
Eqs. pOliep . It reads, 

^'1 (kl, k2) = f C^^i ') - l)-' / ^'^^ ■ ■ ■ ^''^^ 

X •^l6c^^^(ki,k2,qi,-qi, . . . , -Qp; s) Po(gi) ■ ■ •^o('7p) 

(21) 

where J^l^^^^^ = ^abcd^l.d^Ud, ■■■Ud^, and we as sumed 
Gaussian initial conditions together with Eqs. (|8ll2p . 
The factor 2p^) number of possibilities for choos- 

ing 2p points out of 2p -I- 2 and {2p — 1)!! is the number 
of ways 2p points can be connected together by pairs. 

On the other hand, with the relation (jTU]) one can 
also get a formal expression for the cross-correlation 
(*a(k)06(-ki)0c(-k2)). Given the fact that = 
0, such an expression involves the same type of contrac- 
tions of the (paiO^i) factors leading to, 

(*a(k)(^f,(-ki)0,(-k2)) = 2fc(k-ki-k2) 

X Pdb{kl)Pec{k2) (22) 



n"'(k,Pi p„)= 




FIG. 2: Representation of the first two terms of the multi- 
point propagator F'"' in a perturbative expansion. F'"' rep- 
resents the average value of the emerging nonlinear mode k 
given n initial modes in the linear regime. Here we show the 
first two contributions; tree-level and one-loop. Note that 
each object represents a collection of (topologically) different 
diagrams. 




FIG. 3: The one-loop contribution to Gabik, S2, si). The ® 
represents a primordial power spectrum Poiq) with the cor- 
responding "loop" momentum q integrated over with weight 
(27r)^'^ J d'^q . See Q for an explicit calculation of this dia- 
gram. 

with Pab{k) defined in Eq. (11211 . This relation is very use- 
ful since it provides a simple way to compute the three- 
point propagator in numerical simulations out of cross- 
correlations between final and initial conditions (see Sec- 
tion |VT|- 

As expected, the relations in Eqs. ()2m22|) can be 
straightforwardly extended to any order. We thus have 
for the generic case, 

{^ai^)^cA-^l)---Ki-K))c = p!'5D(k-ki...p) 

X ^aL.b, Pb,cAkl)...Pb,cAkp) (23) 

where ( . . . )^ stands for the connected part of the ensem- 
ble average. 



IV. THE LARGE-fc BEHAVIOR OF 
MULTI-POINT PROPAGATORS 

A. The Large-A; Limit of the Two-Point Propagator 

As discussed in the previous section, the two-point 
propagator Gab generalizes gab beyond linear theory and 
thus reflects a key property of the evolved fields. The 
general properties of Gab have been explored in detail 
in Q, but we briefly recall them here to motivate their 
generalization to multi-point propagators. 

Following Eqs. ((gTTU)) . and the definition in Eq. IT^ . 
one can expand the function Gab with respect to the am- 
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0.2 0,4 0.6 



[(Mpc h-i)^] 

FIG. 4: The large-fc limit of the two-point density propagator 
r'^'. Symbols correspond to measurements in numerical sim- 
ulations at redshifts z — 1, 0.5 and z = (top to bottom), see 
text for details. The solid lines correspond to the large-fc limit 
expression given in Eq. ((25}. The linear relation obtained by 
plotting log G vs. fc^ makes evident that the suppression of G 
is indeed Gaussian in the high-fc limit. Moreover, the slope is 
very well predicted by Eqs. (|25I26|I . 



plitude of initial fluctuations, 

Gab{k, Sf, Si) = gab{sf - Si) + G"™P(A:, Sf,Si) + ... (24) 

where G^^°^{k,Sf,Si) is the first nonUnear correction 
term, describing the transition into the nonlinear regime. 
Graphically, this term corresponds to a "one-loop" dia- 
gram (i.e. an integral over Pq), which is shown in Fig. [31 
As nonlinear effects become important Gat is expected 
to decay to zero since they erase the one-to-one corre- 
spondence of modes valid in the linear regime. This in- 
troduces a characteristic scale that describes the decay 
length of the two-point propagator. It was shown in Q 
that this decay can be computed exactly in the high-fc 
limit, where a subset of diagrams is expected to provide 
the dominant contribution. Following a line of calcula- 
tion that we will use again shortly, it was shown that in 
the large-fc limit. 



Gab{k,Sf,Si) = exp [^-Y'^li^"' - e''')^j 3ah(s/ - Si), 

(25) 

where the characteristic decay length is determined by 
the rms velocity fluctuations 

"'^^Jo ^^^^ 




FIG. 5: Example of diagrams contributing to Gab(k) (top) 
and r^^(j*^^(k, ki, k2, ka) (bottom). The dominant contribu- 
tion after resumming all possible configurations is expected 
to come from those diagrams where all loops are directly con- 
nected to the principal line (top) or principal tree (bottom). 
The principal line and tree are drawn with a thick solid line. 
The dominant loops are those drawn by dashed lines, while 
the sub-dominant loops are those in dotted lines. 

In it is shown how to match this result valid for 
kcTy ^ 1 to the low-fc behavior described by Eq. ([24]) . to 
obtain a prescription for its full time and fc dependence. 
This prescription was found to be in good agreement with 
numerical simulations at all scales and different redshifts 
for density and velocity divergence propagators. 

Here, we concentrate on the large-fc behavior of the 
density propagator from growing-mode initial conditions, 

r*^^-' = V^^^Ub — Gil + Gi2 (we will henceforth use both 
G and F^^^ to refer to the two-point propagator). We use 
the algorithm presented in [4] to measure F'-'^' based on 
the cross-correlation property in Eq. . We defer a de- 
scription of the simulations used here until section lVll be- 
low. Figure [4] shows F^^^ normalized by the linear growth 
factor t[11^ — gii + gi2, with gab the linear propagator 
defined in Eq. ([7]); the unusual notation for the growth 
factor is used here to emphasize that it is given by the 
tree-contributions to the two-point propagator, this will 
have a natural generalization for multi-point propaga- 
tors. The figure shows logF^^' vs. logfc^ to emphasize 
the Gaussian decay very well predicted by Eq. (|25)) at all 
redshifts with characteristic scale given by Eq. (pS)) . 

In the following sections we investigate similar proper- 
ties, regimes and measurements to those carried out with 
Gab, to the case of the three-point propagator F^^^ and, 
when possible, extend this to F*^"). 



B. Dominant diagrams and principal trees 

To study the high-fc regime of the propagators, the first 
step is to identify the set of diagrams that is expected 
to dominate the large-fc behavior. This identification is 
closely linked to the concept of principal line and its gen- 
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FIG. 6: Dominant loop contributions for the three-point prop- 
agator r^^'^(ki, k2, ks). The final expression is obtained by 
the sum of an infinite number of such loops and over all pos- 
sible interaction times. In the large-fc limit these terms can 
be computed exactly. 



eralization, the principal tree, that all the diagrams for 
propagators contain. These concepts are illustrated in 
the diagrams appearing in Fig. [5] For instance, in the 
top diagram there is a unique way to go from initial time 
(symbolized by the vertical dotted line) to the final time 
without crossing a (E", thus moving always in the direction 
of increasing time. It's easy to check that this is true for 
any diagram contributing to Gab- This is what we call 
its principal line. 

Similarly there exists a unique tree, the principal tree, 
that joins the initial time to the final time in the bottom 
diagram of Fig. [5l The number of branches of this tree 
is equal to the number of points minus one for the prop- 
agator under consideration (i.e. the number of function 
derivatives in the definition in Eq. (fT7|) ). Again, it's easy 
to check that the existence of a principal tree is general 
and does not depend on the number of loops that are at- 
tached to it, although at a given order different diagrams 
may have different principal trees (see bottom diagrams 
of Fig. [T]) . Principal trees appear then as simple general- 
ization of the concept of principal line initially introduced 
in g. 

We now note that in the expansion series, all diagrams 
can be seen as a principal tree with loops attached to its 
branches. Following arguments developed in Q , we claim 
that the general behavior of the multi-point propagators 
are dominated in the large-fc limit by those diagrams in 
which every loop is directly connected to the principal 
tree. 

This idea follows from the fact that multi-point prop- 
agators are a measure of the correlation between final 
and initial configurations, see Eqs. (|14l22l23p . And dia- 
grams that contain only such loops are more effective in 
the process of correlating the final field with the initial 
conditions (that is, the modes running along the principal 
tree always interact with growing modes coming from the 
initial conditions). In particular, this argument remains 
valid at large-fc where these diagrams can be computed 
exactly, as we will show in subsequent sections. 

As explicit examples, we show in the bottom diagram 
of Fig. [5] in dashed lines the kind of loops that con- 



tribute dominantly to the large-fc limit of F^^^ (while the 
principal tree is represented by solid lines and the sub- 
dominant loops in dotted lines) . In Fig. [HI we show some 
of the dominant loop contributions to F^^^ as they are 
the subject of next section. 

These examples can be used to build a simple picture 
of the different loop contributions. Let us take the case 
of the two-point propagator. Its one-loop contribution 
is shown in Fig. [31 As shown in 4], this contribution is 
dominated by loop momenta of order q ~ cr^^, thus is 
in the high-fc limit {kay 3> 1), most of the momentum 
is fiowing along the principal line of the diagram. Now 
consider the two and higher loop contributions: the dom- 
inant loops, when resummed, give Eq. (PS)) . again saying 
that all dominant loops have momenta of order q ~ 
fiowing in them; this is a consequence of the form of the 
vertices and the shape of the CDM spectrum. 

As far as the subdominant loops are concerned, there 
are two cases at the two-loop level: with one leg in the 
principal line (e.g. as the dotted loop in the top half 
of the Gab diagram in Fig. [5]), and with no legs in the 
principal line (e.g. as the dotted loop in the bottom 
half of the Gat diagram in Fig. [5]). In fact, there are 
two types of diagrams of the latter type, depending on 
whether the two legs of the loop land are on the same side 
of (8) (as in Fig. [5]), or not. The first case corresponds to 
renormalizing the internal propagator, the second to a 
renormalization of the initial power spectrum. However, 
since both renormalizations must be evaluated at large- 
scales, i.e. 9 < cr~^ , their effect on the final result must be 
subdominant. The subdominant loops with one leg in the 
principal line corresponds to renormalizing the vertex, 
and the level of agreement seen in Fig. |4] suggests that 
such vertex renormalizations are small. This corresponds 
to studying the three-point propagator in the limit where 
one incoming momentum is "soft" (of order q < a~^) 
while the other is large. While we are not able to resum 
this case, we shall see that when the whole principal tree 
is in the high-fc limit, the behavior of the three-point 
propagator can be established. 



C. The Large-fc Limit of the Three-Point 
Propagator 

The aim of this section is to obtain the large-fc limit of 
the three-point propagator F^^' defined in Eq. p6| . To 
this end we need to sum over the whole set of dominant 
loop diagrams, 

p(2) _ V- p(2) .271 
^ abc ~ ^ abc, n-loops' ' ) 

n>0 

that is, only consider those diagrams with loops directly 
attached to the principal tree, see Fig. [6] Note that to do 
the summation in Eq. (j27p one needs to take into account 
the symmetry factors of each contribution [2^. Let us 
start by considering the one-loop case in detail first, since 
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it will allow an easy extrapolation when accounting for 
arbitrary number of loops in the large- A: limit. 

Diagrams are computed according to the rules given in 
sectionini but let us be a bit more precise in the notation. 
In Fig. [6l each of the contributing loops carries a wave 
vector q such that q ~ <C fci and similarly g <C fc2. 
It is connected to two segments (not necessarily distinct) 
taken in the set of three - denoted (1), (2) and (3) in 



J 



Fig. [6] - each carrying respectively the wave vectors ki , 
k2 and ks. In the calculation of those diagrams one has 
to sum over all indices, describing density and velocity 
divergence perturbations (we implicitly assume that such 
indices when repeated are summed over). As an example, 
we write down the corresponding integral when the loop 
is connected to branches 3 and 1, 



'"i6l{i3}(^3,ki,k2;s) = jpQ{q)A^q jdsi J ^s' J ds2gad{s - si) 7rfe/(k3, q, kg - q) e^^u^ gfh{si - s') 

'Xlhij{^3 - q,ki - q, k2)giA;(s' - S2)7fe(,„(ki - q, -q, ki) e^'^ui g„rb{s2) gjds') (28) 



r 



where si, S2 are the interaction times where the loop legs 
attach to the principal tree, and s' is the interaction time 
belonging to the principal tree. Furthermore, we assumed 
growing mode initial conditions, so the initial spectrum 
Pab reduces to UaUbPo, see Eq. (|12p. and gab{s)uh = e^Ua- 
It is crucial to note that, as for the calculation of Gab, 
in the high-fc limit the vertex matrix takes the form '4l|, 



7afcc(ki,q,kj - c\)ub 



2q' 



Sac 



(29) 



when contracted with an incident growing mode (such 
as those coming from the initial conditions). The in- 
teraction then leaves the modes along the principal tree 
lines unchanged and carrying the wave vector kj (since 



ki -I- q « ki) all along, 
large-fc limit to Eq. ([28 



As we now show, applying this 
simplifies calculations greatly. 



Let us be general now and consider a loop that joins 
the segment («) and (j). Again, we denote the inter- 
action times for such a loop as si and S2- The loop 
contribution will eventually be given by the integration 
of these times over the time range of the corresponding 
segments. The integral over the segment (i) can be de- 
fined through the respective characteristic functions Xi, 
with Ti{t) = l2{t) = Q{t)Q{s' — t) for interactions after 
the initial conditions but before the vertex of the prin- 
cipal tree (at time s') and l3(i) — Q{t — s')0(s — t) 
for interactions after s' and before the time of interest s 
(s > s' > 0). 



The one-loop contribution is then given by 



.(2) 
abc,{ij} 



-Sij J ds'gad(s - s')7de/(ki,k2,k3)5e6(s')ff/c(s') J dsili{si) J ds2lj{s2) 
X J d'Q.qMqPoiq) 



2g2 2q2 



(30) 



where Sij is the symmetry factor associated to this contribution. It is 4 if i 7^ j (e.g. 2 for each incoming line), and 
it is only 2 for i = j because the time ordering of si and S2 is ignored in Eq. (|30p [29i] . The integral over the angles 
flq is straightforward and leads to, 



.(2) 
abc,{ij} 



J ds'5iad(s- S')7<ie/(kl,k2,k3)gef,(s')5r/c(s') J dsili{si) y ds22'i(s2) kj.kj (T^ e 



2 ^Sl+S2 



(31) 



where cr^ is defined in Eq. ([26|l . Note that, if we could ignore the last two time integrals, the answer would be 
proportional to the tree-level amplitude (given by the first integral). The time integration over si and S2 can be 
done easily for each loop diagram characterized by {ij} = {11}, {22}, {33}, {12}, {13}, {23}. When we collect all the 
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one-loop contributions we get, 



.(2) 
abc, Hoop 



ij} 



ds'gad{s - s')7de/(ki,k2,k3)geb(s')5/c(s') CTv 



ly {kl + 2ki.k2 + fc^) + 2 (e^ - e^') (e^' - l) (ki.kg + ks.kg) + (e^ - e^')' fc. 



(32) 



which when using that ka = ki + k2 makes the square 
bracket independent of s' and leads to the remarkable 
result: 

p(2) __^|_^p(2) fe^_n2 ^oo^ 

ahc, Hoop ~ 2 afcc, tree ^) ' 

We are now ready to consider diagrams with an arbitrary 
number of loops. As we shall see it amounts to exponenti- 
ate the one-loop result above. There are at least two ways 
to obtain such a result. The most direct method, which 
we follow here, consists in computing simultaneously an 
arbitrary number of loops. The cumbersome part in this 
case is to find out the combinatorial weights for each di- 
agram. Here we follow the one-loop strategy discussed 
above and ignore time ordering at first and correct for 
the overcounting afterwards. In Appendix VK\ we provide 
an alternative derivation of this method that takes into 
account time ordering explicitly. A second method based 
on a different resummation scheme is presented in Ap- 
pendix |BJ 



The lesson we must recall from the one-loop calculation 
is the following. Every loop must be directly connected to 
two "principal" segments (i) and {j) of the principal tree 
{i and j go from 1 to 3, see Fig. (S]). Time integration 
then leads to (e'* 1) if i = 1 or 2, and (e* — ) if 
i = 3, and similarly for j. Momentum integration leads to 
— (T^ ki-kj/4, and in addition there is an overall symmetry 
factor Sij (with = 4 if j 7^ i and 2 otherwise). 

The loop resummation in Eq. (P7)) relies on a proper 
counting of each of the diagrams that contribute in the 
large-fc limit. In turn, each of these diagrams can be 
defined by a set of parameters, pij, where each counts 
the number of lines connecting the segments (i) and {j) 
(i ^ j)- For instance, the diagram in Fig. [S] corresponds 
to pn = 1, P22 = 0, pi2 = 1, P13 = 1, P23 = and 
P33 = 1. Each loop line connects to the principal tree 
at time Si. Dropping the time ordering in Si (within 
each of the principal lines), the contribution from such 
an ensemble of loops is given by. 



.(2) 
abc,{pij} 



X (e^' - 1) 



5l 

4 



2 \ Ei<,- Pi] 



2pii+2p22+2pi2+Pl3+!323 



nfcf"n(k,.k,)P- / ds'gUs - s'hdef{ki,k2M)9eb{s')9fc{s') 



(34) 



In this expression the coefficient s^p..j is the symmetry 
factor associated with each set of Pij . For a given ordered 
topology, each loop introduces a factor 4, so that 

^{P..}=2'^'^^-"--- (35) 

The coefficient A4^p.^y on the other hand represents the 
number of times each given diagram has been computed 
because the time ordering has been dropped. Changing 
the time order (within each principal line) allows to auto- 
matically scan the different topologies (diagrams sharing 
the same set oi pij) but computing them several times. 
The symmetry factor pij in M^p.^y is precisely here to 



correct for this multiplicity factor. It is relatively easy to 
see that its dependence on each pij factorizes, so that 

M{p,,}^l[M{p^,) (36) 

i<j 

where A^(pi,;) = 2P"pii\, and M{pij) = pijl Hi ^ j (there 
is an extra factor of 2 when i = j to compensate for the 
absence of time ordering. See Fig. [7] for illustration of 
these symmetry factor countings). The summation over 
Pij is then straightforward and leads to exponentiations, 
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Pll=3 ....-®.. 



^1 



FIG. 7: This figure illustrates the effect of the time ordering exchanges (thick double arrow lines). 

Through such exchanges, the complete set of diagrams that correspond to a given {pij} can be explored. Successive time 
exchanges can however lead to identical diagrams (e.g. left and right diagrams in each panel). The unordered time integration 
then leads to a multiplicity factor for each diagram. For the top panel (corresponding to pi2 — 3), the same diagram is obtained 
each time two loop lines are exchanged. There are P12! of such possible exchanges. For the bottom panel (corresponding to 
pii = 3), the same diagram can be obtained either by the exchange of loop lines, or by the exchange of the end and final time 
of each loop of them (double arrow dashed line). There are thus pn! 2^^^ of such possible exchanges. 



xexp<; -y 



(e*' - 1)' {kj + 2ki.k2 + kl)+2 (e' - s"') (e"' - l) (ki.kg + ka.kg) + ( 



exp 



'^"^S ( . _ ^^2 \ (2 



(37) 



with T^^l^^^^ defined in Eq. This is a truly re- 

markable result. It shows that the whole effect of loop 
summation is encoded in the value of k^ in exactly the 
same way as for the two-point propagator Gab- 

We now compare this result to measurements in nu- 
merical simulations, which will be described in detail 
in section IVII As done for the two-point propagator 
(see Fig. |4]), we test for the Gaussian decay in the 
high-fc limit by plotting in Fig. [8] equilateral configura- 
tions logr^^\k, fc, k) vs fc^, for which Eq. ([57)1 predicts a 
straight line with known slope. We do so for three differ- 
ent redshifts z = 0, 0.5, 1, finding very good agreement in 
all cases with the predictions of Eq. ([37|) . shown in solid 
lines. This validates our resummation scheme. 

Equation (j37p and its generalization to other multi- 
point propagators have important implications for the 
power spectrum and higher-order statistics, that we dis- 
cuss in section IVlIl We note also that a second, faster, 
method to perform the loop resummation is discussed 
in appendix [B) We now consider the extension of these 
results to arbitrary multi-point propagators. 



D. The Large-fc Limit for Higher-Order 
Multi-Point Propagators 

The structure we found for the three-point propagator 
r*-^) is appealing enough to consider its full generalization 
to propagators of arbitrary number of points. The crucial 
property is the extension of the one-loop relation given 
in Eq. ([55)1 . For higher than three-point propagators, the 
tree-order is given by the sum of several diagrams. When 
loop corrections are taken into account, the principal tree 
of each of these diagrams is naturally the corresponding 
tree-order diagram. For instance, there are two of such 
diagrams for F^'*-' shown in Fig. [TJ 
(%>) 

Let us define T, { , , the tree-order contribution to 

(q) .tree ' 

r(P) from diagram of type (g). The value of such a dia- 
gram is obtained after the integration over the si . . . Sp-i 
time variables corresponding to the interaction points is 
done. Generically it can be written as, 

rS,tree = . . . ^ ds^.^T^"^ (ki , . . . , k,_i, 

,si,...,Sp-i) (38) 
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where T^'^^ki, . . . ,kp_i, si, . . . , Sp_i) denotes the value 
of the diagram (g) before the time integrations are per- 




formed. T depends on the vertices and the hnear prop- 
agators, the arguments of which impose time-ordering 
through gabiv) — fo^' ?? < 0. 



[(Mpc h-i)"] 



FIG. 8: The large-fc limit of the three-point density propa- 

(2) (2) 

gator Tl = Tlf^^UbUc, the only density contraction that can 
be measured for growing mode initial conditions, ut — (1, 1). 
The symbols in the figure correspond to equilateral configu- 
rations at redshifts z = 1,0.5,0 (from top to bottom). We 
have normalized these measurements to its low-fc limit T{ ^oc 
given by Eq. (|20p . The figure clearly shows that the measured 
propagator closely follows the large-fc limit given by Eq. (|37)) 
represented by solid lines, once fJ^'^' decays by ~ e~^ from its 
tree-level value. 



The one-loop diagrams built out of the 
(g)— tree diagram involve the same core function 
T*^'')(ki, . . . , kp_i, si, . . . , Sp_i). A one-loop diagram 
joining the segments (z) and (j) is then given by 



— — ^ / dsi... dsp_i T'^''^(ki, . . . ,kp_i, si, . . . , Sp^i) C{y}(ki, . . . ,kp_i, si, . . . , Sp^i) (39) 
^ Jo Jo 



f (k.-k,)(e^ 



(40) 



where s^.i and Si^2 are respectively the starting and end- 
ing time of segment (i) of the principal tree which carries 
momentum k^. Thus the times Skj belong to the set 
{si, . . . , Sp_i}. This is already the result of integrating 
over the additional two time variables that correspond to 
the two extra interactions due to the loop (see the step 
in going from Eq. ([?T|) to Eq. (|32p . for example). 

From these results it can then be shown that the sum 
over of C{ij} is \ki + ■ ■ ■ + kp_i|2(e" - 1)^. The 

proof can be obtained recursively as a simple successive 
application of the result found for F'^' . This mechanism 
is illustrated in Fig. O It leads to (kp = ki + • • • + kp_i), 

-'"(gj.lloop = ^^^^^ ~ -'-)^-'"(g),troo- (^l) 

As shown in appendix [B] the scheme for resummation 



over all loop terms can be applied once again, leading to 
the generalization of Eq. ([57)1 , 

F(^)=exp[-^(e^-l)2] Ftl (42) 



V. ONE-LOOP RESULTS FOR THE 
THREE-POINT PROPAGATOR 

Having derived expressions for large scales, the tree- 
level expression in Eq. p^ . and small scales (high-fc 
limit), we now discuss what happens in between. This 
can be achieved by computing the next to leading order 
correction in an series expansion in terms of initial fields, 
i.e. p = 1 in Eq. (j2ip . Diagrammatically it corresponds 
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\ \ 

FIG. 9: Sum of the one-loop contributions described by C(ij), defined in Eq. (I40p . The contributions involving the wave 
vectors ki, kj and their sum is singled out. The first three diagrams correspond to the one- loop contribution involv- 
ing an extra wave vector, k' (dashed line). By the high-fe summation rules described before Eq. (|34p . they sum to 
k' ■ [(k, + kj){e'^ - e^i) -f k,(e=i - 1) + kj(e''i - 1)] [e"^ - 1) = k' • (k, + kj)(e"=^ - 1)^ The next three and last three di- 
agrams correspond to "internal" loops (i.e. loops attached to the principal lines with momenta k;, kj and their sum). As 
shown for the computation of r'^-*, these six diagrams sum to |ki -I- ^1^(6"^ — 1)^. By the summation rules, these two results 
can be expressed by two contributions to C(ij) with one less branch, as shown by the two diagrams in the last line. The same 
procedure can be applied recursively until the expression of C(ij) reduces to that coming from a single branch, leading to the 
final expression. 



to calculating the one- loop diagrams from Fig.[6l accord- 
ing to the rules given in section [ill There are six types 
of such diagrams, depending on where each of the loop 
lines is attached to in the principal tree. Depending on 
this, they can be labelled as {33}, {11}, {22}, {12}, {13}, 
and {23}. Each of these cases is depicted in Fig. [51 
In Eq. we wrote down explicitly the one-loop inte- 

(2) 

gral expression for F^^^^ {13} ■ ^ similar procedure can be 
followed to obtain expressions for the rest of the one-loop 
contributions. The difficulty in their actual computation 
resides in the integration over loop momenta. The an- 
gular integrations are independent of Po{k) and can be 
done analytically. Such integrations involve two sets of 
functions, one built out of 



with, 



klkj + q^ + q^{kl-kl-k2^) (45) 



1 



|ki+q|2 

and the other built out of. 



71" log / (fci +qf 



kiq \{ki-qy 



|ki+qP|k2-q|- 
7rlog(>V+/>V-) 



kaqVWW+q^ 



q^{kl-kl-k2^) 



(43) 



(44) 



and fc3 = |ki -|-k2 1 . The explicit expression of the final re- 
sult for arbitrary dependence on configuration iki, /c2, ^3) 

is too long to reproduce here. However, we have checked 
(2) 

that ^jq^p (ki , k2 , g) follows the expected large-A: 
limit. In the \ow-k limit, it is interesting to note that 
after the angular integration ra(,c,iioop(ki, k2, g) behaves 
as l/q^. 
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The most growing contribution (i.e. largest power in growth factor) in this limit is, 

"4901(fc^ + 32879fc]^ + 231478fc^fc| + 32879fc| 



rl'Lop(ki,k2,.) 



2T5 (ki • k2 



18865 226380A;fA;| 
9552(fc? + kl) 2 12409 



18865 fc2fc2 ^ ' ^' 56595 kfk? 



<aHr), (46) 



with given in Eq. (|26p and a(T) the linear growth factor. Note that its angular dependence can be different from 
that of the tree order expression, depending on configuration. 

For the complete result (but still most growing) valid at all scales we restrict to particular triangle configurations. 
For equilateral configurations (ki = ^2 = fcs = fc), 

(5754fc*^ - 33521g2/c'' + AUAlq'^k'^ - 27039g^fc2 + 5175g^) 
-3 {k^ - q^f (2174fc4 + 441392fc2 _ 1725^4^ ^/^^^k^T^log (||:^) } , (47) 

where W^'^"' correspond to the value of W± for an equilateral configuration. For colinear configurations {k^ = 2ki = 
2k2 = fc), 

p(2) coiin _ [ p ( \ / 57171840fcg - 69330880g2fc6 + 49644768g^A:'^ - 13102080g^fc2 + 2349765g^ 
i'"°°P " J \ 3179151360fcV ^ 

31og(w|'VW-'') (g' - 4fc2)^ log(>vi'Vwi'^) (g' - 4fc2)' (2fc4 + 7q2fc2 - 9q^) 
275968fcgVl6A;4 - 7q2fc2 + ^4 29568fcg5y4fc4— ^2^-^^ 

(ll46fcio + 3259g2fc8 - 7539g'*fc6 + 4191g6fc4 - I363q^k^ + 306q^°) log ({1^) 

1892352A;V 

(g2 - 4fc2)^ (976064fc6 - 253104g2fc4 „ ijsUAq^k^ + 19563^^) log ( g^+g^ ^ 

1695547392fcV ' 

where >v£' = W±(2fc, 2k, k) and W^^^ = W±(2fc, A:, 2fc). 

I 



Equations (|47|) and ((48l) . fully specify the most grow- 
ing contributions to the density propagator, their time 
dependence being an overall scaling by a^(r). 

As examples we show in Fig. [10] the dependence 
of T\ ^j^^p(ki, k2, 9) on q for particular triangular con- 

figurations, with Fj^ j^j^^p(ki, k2, g) defined such that 

rS-ioop(ki,k2) = /F(^ji„„p(ki,k2,(?)Po(g)g^dg.^ For 
small values of q, corresponding to the large-fc limit, we 
recover the tree-order result. For large values of g, the 
integrand behaves similarly, but with a different ampli- 
tude. 



cross-correlation property, Eq. (|22|) . allows us to easily 
measure the three-point propagator in numerical simula- 
tions. In fact, this can be straightforwardly extended to 
other multi-point propagators using Eq. (|23p . 

In practice initial conditions in simulations are gener- 
ally set in the linear growing mode in which case Eq. (j22[) 
becomes, 

(^'a(k)^o(-ki)(5o(-k2)) = 2fc(k-ki-k2) 

X ri2)(ki,k2)Fo(fci)Po(fc2), (49) 



VI. THE THREE-POINT PROPAGATOR IN 
NUMERICAL SIMULATIONS 

We now show in more detail comparison of our re- 
sults against measurements in numerical simulations for 
the three-point propagator. As mentioned before, the 



where we used the shorthand notation Fa = ^abc^bUc 
with Ua = (1,1) and </'a(k) = Uai5o(k), and the initial 
power spectrum obeys Po(fc)5D (k-f k') ~ (5o(k)(5o(k')))- 

f2) 

Therefore, Fa can be measured by implementing the 
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FIG. 10: Momentum dependence of the three-point density 
propagator integrand T^^^j^^ (ki,k2,(j') for different config- 
urations. The sohd line denotes equilateral triangles (with 
fci = ^2 = fcs = 1); the short dashed line colinear triangles 
(with fci = fc2 = 0.5, ks — 1); and the long dashed line trian- 
gles with fci = 0.816, k2 = 0.3, — 1 (in all cases units are in 
h/Mpc). The y-axis has been normalized to the large-fc limit 
expression given in Eq. (|33p . Notice how this limit is reached 
for small values of q, as expected. 

following algorithm based on Eq. (|49l) . 

xSai-kM-k,), (50) 

where the sum runs over Fourier modes in the |ki| 
bin, kj in the |k2| bin and k; in a bin jkaj such that 
|ki — k2| < Iksl < ki +k2, and N is the number of terms 
in the triple sum. 

In writing Eq. (|5(I1) we have assumed that the three- 

(2) 

point correlator Fa does not depend upon the orien- 
tation of the triangle but only on the magnitude of its 
sides, as expected from statistical isotropy. Additionally 
one must be aware that the estimator in Eq. (|50p is only 
strictly valid when the initial field 5o is Gaussian (oth- 
erwise, one must go back to the definition Eq. (|16p and 
implement the functional derivative, see for such de- 
tails in the case of the two-point propagator). 

In this paper we focus on the statistics of the density 

(2) 

field, and thus we only present measurements of F . Ob- 

(2) 

taining Fj requires the measurement of the local veloc- 

(2) 

ity divergence, while measuring other components of F^^^ 
requires numerical simulations with a mixture of growing 
and decaying modes in the initial conditions. 

To test our theoretical results from sections IIVI and 



rVl we use Eq. ([50] to measure the three-point propa- 
gator in N-body simulations using a set of 50 realiza- 
tions, with each realization containing Npar = 640'^ par- 
ticles within a cubic volume of side L^ox = 1280/i~^Mpc. 
The total comoving volume of our simulations, approxi- 
mately 105 (h^^Gpc)^ , is large enough to minimize cos- 
mic variance. Cosmological parameters were chosen as 
= 0.27, = 0.73, Qb = 0.046 and h = 0.72. The 
corresponding initial power spectrum had scalar spectral 
index rig = 1 and was normalized to yield erg = 0.9 
when linearly evolved to z = 0. The simulations were 
run using Gadget2 [2^ with initial conditions set at 
Zi = 49 usirig 2nd order Lagrangian Perturbation Theory 
(2LPT) [3, m. The output times were &t z = 1,0.5,0. 
At each output time we measured the two and three-point 
propagators by cross-correlating the density field at the 
desired time with the initial density field at Zi — 49 (as 
described above and in y). 

Figure [TT] shows the density three-point propagator 

(2.) 

r\ ' for equilateral configurations (top left panel for z = 
and bottom panels for z = 0.5, 1, as labeled) and colinear 
configurations (top right panel, z = Q). In each panel the 
symbols with error bars correspond to the measurements 
in the N-body simulations, the long dashed line repre- 
sents the one-loop result from Eqs. (|47l48p and the solid 
line the large-fc limit from Eq. (|37p . In all cases the errors 
displayed are for the mean of the ensemble obtained from 
the scatter among the 50 realizations. For simplicity we 
only show the more familiar configurations, correspond- 
ing to colinear and equilateral triangles, although we have 
measured all configurations. 

The comparison in Fig. [11] shows overall that both 
large-fc and \ow-k regimes are well modeled by the ana- 
lytic results for different triangle configurations and red- 
shifts. This establishes a good basis from which a match- 
ing ansatz can be constructed to describe also interme- 
diate scales, as done in [3| for the two-point propagator. 

Let us now point out some limitations in the measure- 
ments that help understand better the significance of the 
level of agreement seen in Fig. [TT] The asymptotic limit 
to the tree level value at large scales is difficult to measure 
since few configurations of a given shape are available at 
low k, which gives rise to relatively large cosmic variance. 
To partially mitigate this one is tempted to increase the 
bin size in Fourier space. However, a larger bin size in- 
troduces a mixing of triangles of slightly different shape, 

(2) 

which have different T\ . For equilateral triangles, in- 
creasing the bin size introduces slightly off equilateral 
triangles which bias the result to higher values as we show 
below, while the opposite is true for colinear configura- 
tions. Therefore a trade-off takes place between reducing 
the bin size to obtain unbiased results but at the expense 
of increasing the cosmic variance. We choose a bin size 
5k — 0.005 ft,Mpc~^ ~ fcp (the fundamental mode of the 
simulation volume) for k < 0.15/iMpc~^, 2kY that up 
to < 0.4/iMpc~^ and 4fcF for smaller scales where the 
impact from offset configurations is negligible. This is 
the reason why the error bars diminish by a factor of 
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a few at e.g. k ~ 0.15/i~^Mpc for equilateral triangles. 
For colinear triangles, the change in binning happens at 
slightly different scales, evident from looking at the top 
right panel in Fig. [TTl 

The bias of the r(2) estimator can be understood and 
quantified using the tree level results as follows. Since 
ka = ki + k2, equilateral triangles have an angle of 120° 
between ki and k2. Configurations slightly off this will 
give a tree level result, from Eq. (|20|) . equal to 

r(^Le(fc, k, cos(120° + <5)) « y - ^<5, (51) 

where 2/7 is the "exactly" equilateral value. The off- 
equilateral configurations with larger have more sta- 
tistical weight (more Fourier modes), and they are rep- 
resented by negative S. Therefore a bias towards higher 

(2) 

values of Ti is introduced by increasing the bin size. 
Conversely, the slightly off-colinear configurations would 
bias the measurement below the exactly colinear case be- 
cause 



Finally, it is easy to show that the size of error bars 
on ri^-* scales as -D+^, as for the bispectrum [2l|. This is 
nicely recovered in Fig. 1111 when one compares z = 1,0.5 
and 0. 



VII. APPLICATIONS TO THE POWER 
SPECTRUM AND BISPECTRUM 

A. Reconstructing the Power Spectrum from 
Multi-Point Propagators 

We now discuss the relationship between multi-point 
propagators and nonlinear corrections to the power spec- 
trum, showing that the power spectrum Pat can be recon- 
structed by gluing together T^p^ contributions, as briefly 
mentioned in Eq. (|18p above. 

A formal expression for Pab{k) can be written using 
the expansion in Eq. ^ as. 



r(%,,(fc, k, cos(0° + 5))^2~ ^<5^ (52) 

resulting in a negative estimator bias. 

The empty circles in Fig.llllcorrespond to the measure- 
ments "corrected" for the estimator bias assuming that 

(2) 

the configuration dependence oiT\ , slightly away from 
the exactly colinear or equilateral configurations, is given 
by the tree-level expressions above (with 5 kp /k) . For 
equilateral configurations the agreement with the one- 
loop prediction at intermediate scales is improved, par- 
ticularly at z = 0.5 and 1. For colinear configurations 
the correction is negligible. 

These estimates assume tree-level perturbation theory, 
when loop corrections are included the situation is ex- 
pected to changed somewhat. As it is well known in the 
bispectrum case [2l[, at low redshift the dependence on 
configuration is suppressed, so the estimator bias is ex- 
pected to be larger at higher redshift where mixing dif- 
ferent triangles leads to more bias. This helps explain 
why at higher redshift (bottom right panel in Fig. [TT|) . 
the intermediate scale measurements done with bin size 
2A:f are more above the predictions than at low redshift. 



SD{ki+k2)Pab{kl,s) EE (*„(ki,s)«'fc(k2,s)) 

= J2 (*i"^nki,.)*^^(k2,s)) 

(53) 

Then, for a given choice of indices ni and n2 one has to 
compute the ensemble average of ni + n2 factors 0(qi), 
following the field expansion in Eq. (fTU]) . 

Since we assume Gaussian initial conditions, this joint 
ensemble average is then a sum of various terms, each 
of them the product of two-point correlators (Wick's 
theorem). Each of these terms can be labelled with 
indices r, s, t, where r is the number of connected pairs 
within the first ni fields (q^), s is the number of 
connected pairs within the last n2 fields (q^, and 
t is the number of "mixed" pairs connecting fields 
from the first rii and the last n2. Obviously one 
has ni — 2r + t and n2 — 2s + t. Let us define 

{(j)c, (qi) . . . (i)c2r+M2r+t)■Ad^ (q'l) • • ■ <l^d2.+t{<^2s+t)) r,s,t 

as the subset of the terms appearing in 

{(j)c^ (qi) . . . (j)c2r+M2T+t)(j)d, (q'l) • ■ ■ (/-d^^+t (qL+t))^ 
that correspond to a given r, s, t triplet. 



Then the discrete sum in ()53p can be written as a sum over r, s and t as. 



(*a(ki)*f,(k2)) = XI / 'i^^i • ■ ■ d^q2r+t d^q'i . . . d^q^.+t ^D(ki - qi...2r+t) fc(k2 ^ q'i...2.+t) 

'^iX^%r+t (qi' • ■ ■ ' q2r+t; »7)-?^MiT.dL+t (q'l' ■ • • ' qL+t; '7)(<^ci (qi) ■ • ■ 0c2.+t i<=i2r+t); (j^d^ (q'l) ■ • • <t)d2.+t (qL+i))^ , t 
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FIG. 11: Measurements of the three-point density propagator T{ against analytical predictions for different triangle configu- 
rations and redshifts. Top Left Panel: equilateral configurations at 2: = 0. Top Right Panel: colinear configurations at z — 0. 
Bottom Panels: equilateral configurations at z — 0.5 (Left) and z — 1 (Right). In each panel the long dashed line corresponds 
to the one-loop prediction given in Eqs. (|47I48|) while the solid line is the high-fc limit derived in Eq. (|37|l . At large scales the 
measurements reach the tree-level result well known from standard perturbation theory (horizontal short dashed line), given 
by Eq. (|20|l . The empty circles in the panels of equilateral configurations are the measurements after the correction due to the 
estimator bias as described in the text. In addition, the sudden change in the size of error bars as a function of scale is due to 
the change in binning, see text for details. 



where we have omitted the time variable rj for clarity. This expression can be further written as, 

(vl/,(ki)*.(k2))=^^^t!(2r-l)!!(2.-l)!!(^2" + *) (^^' + *^ J d'V • • ■ dV dVi • • ■ dV. dV/ ■ • • d^qj' 

x^r+*) (q'Z, . . . , q;', qi, -qi . . . , q., -q.; v) ^t^'^ (-q'l', ■ • • , -q", q'l, -q'l • • • , q^, -q^; v) 

xPo(<Zl) . . . Po(gr) Po{q[) ■ ■ ■ Poiq's) Poiq'l) ■ ■ ■ Po{q't') (55) 
where ^a"^ = ^aci...c„'Uci . . . Wc„ and Pq is the initial power spectrum from Eqs. (j8ll2p . Remarkably, one can recognize 
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FIG. 12: Reconstruction of the power spectrum from multi- 
point propagators, where a large shaded circle represents all 
possible loops that enter into the fully nonlinear multi-point 
propagator (see Fig. [2] for its perturbative expansion) . The 
crossed circles represent initial power spectra. The sum runs 
over the number of internal connecting lines, i.e. the number 
of such circles. It is important to note that each term in this 
sum is positive. 



in this expression a sum of products of r*-*-* functions, see Eq. pT|) . leading to the formal expression, 

Pafc(k, V)=Y.^-J d'qi • • ■ d'qt <5D(k - qi...t) ri*' (qi, . . . , q*; 77) (qi, . . . , q^; 77) Po(gi) • • ■ Pofe), (56) 



where we have introduced the shorthand notation, 

ri*)(qi, • • • , qt) = r(*2,...c,(qi, . . . ,cit)uc, . . .u,^, (57) 
and used the following property, 

^a.bi...b„ (''I' ■ ■ ■ ' ^n) — ^abi...b„ (" ^i, . . . , — k„) 

Thus, we have shown that the power spectrum can be 
written as a summation over contracted n— point prop- 
agators. This sum is diagrammatically represented in 
Fig. [121 and corresponds to a contraction of F^"^ di- 
agrams when the incoming lines are glued together to 
form initial power spectra. The combinatorial factor t\ 
in front of Eq. (j56|) is the number of ways of contracting 
two {t + l)-point propagators. 

It is interesting to note that for a = b, the result is 
a sum of positive terms. For example, for the density 
spectrum (a = b = 1) up to one- loop corrections we have, 

P{k) = [F(i)(fc)]2Po(fc) 

+ 2 J A[r('Hk-q,q)]2Po(|k-q|)Po(g). 

(58) 

This resummation scheme is then very attractive in 
practice, because, unlike standard PT, the subsequent 
contributions add constructively and correspond to a sum 
over "bumps" peaked at increasingly higher fc's. We en- 
counter here a similar situation as the one presented in 
the original formulation of RPT [1, @, although the re- 
summation scheme differs in detail. For example, here 
it is important to note that we are led to consider an 



infinite number of objects, the n— point propagators, a 
property at variance with what a renormalization group 
approach suggests [33|. This is unlike in 0, Q where 
one deals with a few objects, e.g. the two-point propa- 
gator and the vertex (in this regard, note that Eq. ((56)) 
automatically incorporates vertex renormalization). We 
have seen though that the multi-point propagators share 
common, and useful, properties. In particular, their ex- 
pression for the high-fc limit (fccr„ 3> 1) provide a useful 
approximation to compute multi-loop corrections to the 
power spectrum, i.e. the mode-coupling power spectrum 
in Eq. (HID. 

We finalize this section by noting that the relation in 
Eq. (|56p can be also obtained by considering an expansion 
of the field in terms of generalized Wiener-Hermite 
functionals of the initial fields (pa, see The multi- 

point propagators defined in Eq. (flT)) would then cor- 
respond to the kernel coefficients of this expansion. In 
[r2j it is also introduced the formal expansion to address 
the more general case of non Gaussian initial conditions. 
Note however that this approach does not provide a pre- 
scription to compute these kernels (e.g. a resummation 
scheme such as the one derived in Sec lIVI or in [5| for 
F(2)). 



B. Reconstructing the Bispectrum from 
Multi-Point Propagators 

We now extend these results to the bispectrum. Gen- 
eralizing Eq. (|56p it is easy to show that the bispectrum 
can be formally written as. 
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FIG. 13: Reconstruction of the bispectrum from multi-point 
propagators. The crossed circles represent initial power spec- 
tra. The sum in Eq. (|59|l runs over the number of connecting 
lines between each of the emerging modes, e.g. that cross 
each of the dashed half lines. 



(*a(kl)^'f,(k2)«'c(k2)) = J2 



r + s 
r 



s + t\ ft + r 



'■IsW. J d^qi . . . d=^q^ d^q; . . . d^q^ d^q'/ . . . d 



s / \ t 

xfc(ki - qi...r - qi...J (5D(k2 + qi...^ - qi...J SvCks + q" j + qi...^) 

xr(:+^) (qi, . . . , q,, q'l, . . . , q^) ri^+*' (~q'i, . . . , ~q'„ q'/, . . . , qj') 

X r^--) (-q'Z, . . . , -q;', -qi, • • ■ , -q.) Po(gi) • • ■ Po{qr) Po{q[) ■ ■ ■ Po{q's) P^{qT, 



Poiq't')- 
(59) 



This sum is diagrammatically represented in Fig. [T31 We see that it runs over the number of lines that connect each 
side of the diagram (with the constraint that at most one of the indices r, s or t is zero, otherwise we would have 
a disconnected diagram). The leading order (tree) contribution is then obtained forr = s = l,i=:0 (plus cyclic 
permutations), up to one-loop corrections (in square brackets) we have 



B{h,k2,h) = 2r(2)(ki,k2)r(i)(fci)r(i)(fc2)Po(fci)i^o(fc2)+cyc. 

+ [sj d3qr(2)(ki-q,q)r(2)(k2+q,-q)r(2)(q-ki,-k2-q)Po(|ki-q|)Fo(|k2+q|)Po(g) 

+ 6 y d3qr(3)(_k3, -k2 -I- q, -q)r(2)(k2 - q, q)r(l) (k3)Po(|k2 - q|)Po(<z)Po(fc3) + cyc 

I 



(60) 



Note that having resummed the multi-point propagators 
means that many of the one-loop corrections in standard 
PT are already encoded in F^^'^ and thus the number of 
one-loop diagrams is reduced. For the power spectrum 
we have one instead of two diagrams, for the bispectrum 
we have two instead of the four in standard PT [llj . 

It is useful to compare the structure of Eqs. ((58|) 
and ([50]) . We see that the one- loop corrections to the 
power spectrum depend on the initial power spectrum Pq 
through a convolution with the three-point propagator 
F^^), which determines the large-scale (tree-level) bispec- 
trum. The two-loop correction to the power spectrum in- 
volves a similar convolution with F^^^ , which determines 
the large-scale trispcctrum, and contributes to the one- 
loop bispectrum. This pattern continues to higher or- 
ders, demonstrating that in order to extract the most 



information about the initial power spectrum Pq, it is 
advantageous to simultaneously measure the power spec- 
trum and higher-order spectra at large scales and include 
these relationships when doing cosmological parameter 
estimation. 

As a preliminary application of these results, we com- 
pute the reduced bispectrum Q defined by 



P(fci, k2, fcs) 



P(fcl)P(fc2) + P(fc2)P(fc3) + P{ki)P{ki 



(61) 



where we use one-loop results for both the power spec- 
trum and bispectrum from Eq. ([58]) and Eq. ((60)) . respec- 
tively. Since we don't yet have a full prescription for the 
multi-point propagators valid at all scales, we use their 
high-fc limit expressions, Eq. (H^ modified as follows. 
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i - (1) . ^ tree, [^Z) 
tree V^py 

which reduces to Eq. (j42p in the high-fc hmit, and approx- 
imately incorporates (through the fuU scale-dependence 
of the two-point propagator T^^'> found in the fact 
that at low-A: propagators decay slower than their high- 
fc limit indicates, as shown in Fig. [11] by comparing the 
one-loop to the high-fc expressions at low-fc. 

Figure [T3] shows the result of this calculation (solid 
lines) compared to measurements of the bispectrum in 
our numerical simulations (symbols with error bars) . Be- 
cause we only compute up to one-loop corrections, the 
range of scales available to test our predictions is lim- 
ited. At low-fc, the deviations from tree-level FT are 
small, and error bars increase. At high-fc two-loop cor- 
rections become important (see e.g. Fig. 3 in [6] for the 
case of the power spectrum). We concentrate on z = 
where the deviations from tree-level are largest and the 
application of standard FT is least successful. The scales 
shown in Fig. [T3] are those where the deviations from 
tree-level are largest while we seem to be still within the 
validity of the one-loop approximation. We see that the 
level of agreement with simulations is very encouraging 
despite the approximate nature of Eq. , correctly de- 
scribing the loss of dependence of Q on triangle shape. 
More work is needed to explore the robustness of these 
predictions, but it is clear that the results found in this 
paper on multi-point propagators will be useful for a vari- 
ety of statistics. A detailed discussion of the bispectrum 
will be presented elsewhere. 



VIII. CONCLUSIONS 

We generalized the notion of two-point propagator to 
arbitrary multi-point propagators. We derived the one- 
loop corrections to the three-point propagator, and were 
able to resum the dominant class of diagrams contribut- 
ing to multi-point propagators in the large-fc limit, ex- 
tending what was derived before for the two-point prop- 
agator in Q. Our results show that multi-point propa- 
gators decay with the same rate as the two-point prop- 
agator in the large-fc limit and are proportional to their 
tree-level values. We have extended the algorithm pre- 
sented in Q to measure the three-point propagator in 
numerical simulations, and showed that our high-fc limit 
resummation agrees well with both two and three-point 
propagators at nonlinear scales. 

Furthermore, we showed that the multi-point propa- 
gators can be considered building blocks out of which 
one can construct arbitrary n— point spectra. In addi- 
tion, nonlinear corrections to the power spectrum depend 
on the initial spectrum through convolutions with multi- 
point propagators, which in turn determine the large- 
scale (tree-level) polyspectra. This gives a direct physical 
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FIG. 14: The reduced bispectrum Q at z = in tree-level 
PT (dashed line), one-loop RPT (solid line) and measured in 
numerical simulations (symbols with error bars) as a function 
of the angle between ki and k2. 

connection between nonlinear corrections to the power 
spectrum at small scales and higher-order correlations at 
large scales. 

We used our results on the resummation of multi-point 
propagators to reconstruct the power spectrum and bis- 
pectrum up to one-loop level, and showed that the pre- 
dictions for the reduced bispectrum agree well with mea- 
surements in numerical simulations. This is encouraging 
news for the calculation of higher-order correlation func- 
tions in renormalized perturbation theories. 
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APPENDIX A: TIME-ORDERING AND 
EQUATION ([37]) 

Here we provide another derivation of the resummation 
of the three-point propagator that exphcitly takes into 
account time-ordering. An n-loop diagram is specified by 
the number of vertices nt (i = 1, 2, 3) along the principal 
lines with momenta ki, i.e. 2n = ni + n2 + n^, ni = 

2pil + Pl2 + Pl3, "2 = 2p22 + Pl2 + P23, ri s = 2 p33 + 

Pi3 + P23 , where pij are defined as in section IIV CI (note 
that n = ^i<:j Pij)- Here we start from scratch, i.e. 
to construct such an n-loop diagram we start from the 
principal tree, to which we attach branches to principal 
lines of momentum ki each with its initial condition (as 
in Fig. [1]) and we proceed to take expectation values. 

Consider one of the principal lines at a time, say (3) 
carrying momentum ^3. There are 713 time-ordered ver- 
tices along this line, and its symmetry factor is 2"=* be- 
cause branching is asymmetric (which leads to Eq. (j35p 
when all three principal lines are included) . Out of these 
71-3 branches one must select 2^33 for averaging among 
themselves (thus making ^33 loops), out of the remain- 
ing pi3 -|- P23 one must select P13 for averaging with pi^ 
branches attached to the principal line (1), thus creat- 
ing pi3 loops, and the remaining P23 average with P23 
branches attached to (2). The number of ways of doing 
this is. 



n3 

2P33 



/ Pi3 + P23 

\ Pl3 



713! 



(2p33)!pi3!p23! 



(Al) 



and the analogous combinatorial factors also apply to the 
other 2 principal lines. The process of taking expectation 
values introduces the following factors: {2pii — 1)!! is the 



number of ways of averaging pairs (Gaussian fields) out 
of 2pn initial conditions, while there are Pij \ ways of con- 
necting branches {i ^ j) coming out of principal lines 
{{) and {]). This leads to the combinatorial factor, 



(2pn - 1)!! (2p22 - 1)!! (2^33 - I)!!pi2!pi3!p23! (A2) 

Putting together this with Eq. (|Aip and its two permuta- 
tions for the two other principal lines (and the symmetry 
factors for each line) we have a total overall factor of 



2ni+n2+«3 773! 



(A3) 



Finally the time-ordered integrals give, e.g. for principal 
line (3), 



dsic"! / ds2e'* 



ds„3e^"3 



(A4) 

where s is the final time and the principal tree splits 
at time s'. For principal lines (1) and (2) one obtains 
analogous factors but with the corresponding time de- 
pendences, i.e. (e'' — 1). Therefore, we see that the 71,^! 
factors in Eq. (jA3[) get cancelled by time-ordering (as ex- 
pected from the discussion in section flV C) see in partic- 
ular Fig. [7|). Then, after momentum integration, which 
gives —a'^kikj / A for loops attached to principal lines (i) 
and (j), we can write this equation taking into account 
time ordering as follows. 



= n ^''-5'^'" n^^^^^T^ /'d.',..(.-.')7.e/(k„k2,k3),..(.')5/e(.') 

X f - Ij ( - eM , (A5) 



which when summed over the pij, and using that k3 = 
ki-|-k2, gives the tree-level value times a Gaussian factor, 
i.e. Eq. (1271). 



APPENDIX B: RESUMMATION IN FOURIER 
SPACE SHELLS 

The results presented in the main part of this paper 
rely on a direct resummation of multiloop contributions. 
In this appendix we present another method that makes 
use of the mathematical structure of the loop correction 
terms in the large-fc limit. The idea is to partition Fourier 



space into finite volume elements. Although there is no 
need to actually specify the partition we will be using, one 
simple way of constructing an explicit one is to decom- 
pose the volume into concentric shells that we will denote 
i, j, . . . , with the i-th shell corresponding to q G [Ai_i, A^] 
and being of width ^A^. That means that the integration 
over Fourier modes, / dq'^, becomes 47r^- AfSAi assum- 
ing that the integrand depends on the norm of q only. Let 
us now explore the loop contributions for the two-point 
propagator Gab m this framework. 

According to the results presented in section IIV[ see 
e.g. Eq. (PO)) . the contribution to the one-loop correction 
coming from shell i to Gab{k), 5^Gab{k), reads in the 
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large- fc limit 

where we used the following definitions, 



(Bl) 



£(s2,si,fc,A) - -M!e-ie^=Po(A) (B2) 

C{k,K) = / dsi / As2 C(s2,si,k,K) 
Jo Jo 

= _M!(e^_i)2p„(A). (B3) 

Multiloop corrective terms exhibit then a truly remark- 
able property. Let us start with a two-loop contribution. 
In all of these calculations the linear propagator factors 
along the principal line can be factorized in the large-fc 
limit, due to the property that gabis — s')gi,c{s') — gads)- 
Thus, we simply have for the two-loop contribution to 

Gab, 5^'^Gab-, 

5'''Gab{k) = gab C(k,Ki)5K, C{k,Aj)SA,. (B4) 

This factorization property can be generalized to any or- 
der. As a result, for a collection of loops whose indices 
are in set X, the multi-loop correction term, 5^ Gab, reads 



5^Gab{k)^gab\{[C{k,K,)5K, 



(B5) 



We can then resum all these contributions, e.g. resum 
S^Gab for relevant sets I of shell indices. Let us denote 
Gab,Ai (k) the resummed value of the propagator when 
all contributions of shells up to i have been taken into 
account, that is, assuming there is a cutoff for wavemode 
q at g = Aj. 

The result is formally given by a sum of terms of the 
form S^Gab{k) where the sum runs over all possible sub- 
sets X of {1, ... , i}, e.g., 

Gab,A. = S^Gabik), (B6) 

IC{1, ...,»} 

and when the width of each shell is made infinitely small. 
The subsets of {0, ■ ■ ■ ,i} can be split into the subsets of 



{0, . . . , z — 1} and the subsets of {0, ...,?} that contains 
j. As a result this sum can be rewritten. 



ICl,...,i 



S^Gabik) = Y ^^Gab{k) 

ICl,...,i-l 

+ S^^'^'^Gabik). (B7) 

ICl,...,i-l 

Because of the factorization property, (jB5p . this decom- 
position then reads 

Gab,Ai ~ Gab,A,-i + Gq{,,A;_i £(fc, Aj)JAi . (B8) 

In the continuous limit it leads to the equation, 

SGabik, A) = Gabik, A) C{k, A)dA. (B9) 
The solution of this ODE is Eq. ((25|) for the initial con- 
dition Gab{k, 0) = gabik). 

This procedure can be simply extended to multi-point 

(2) (2) 

propagators. Let us denote r)jj^ ^. the expression of r)jj^ 
when modes above Ai are cut out. In a construction 
similar to that of the two-point propagator, the structure 
of the loop correction is such that. 



^M'c(kl,k2,k3) = r:,tc,tree(kl,k2,k; 

iei 

It them implies that 

^abcA. (kl, k2, k3) - r^^l^^_^ (ki, k2, ka) = 

ri'U(ki,k2,k3)/:(fc3,A,)M,. (Bii) 

In the continuous limit it leads to a simple ODE that can 
straightforwardly integrated into (|37p taking into account 
the initial condition, 



(2) 



(BIO) 



p(2) ^ p(2) 

abc.A—0 abc, tree' 



(B12) 



This construction extends to any multi-point propaga- 
tor, or more precisely the procedure extends to the com- 
putation of the loop corrections that apply to any distinct 
tree structure of tree order multi-point propagators. 
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